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Optimum ranges of the inertia ratios, the spring constants, and the 
damping constants have been obtained for the design of a gravitationally 
oriented two-body satellite with satisfactory over-all damping performance. 
In the case of viscous damping, optimum damping constants can be simply 
chosen from diagrams of complex root loci. It is found impossible to con- 
vert the optimum viscous damping constants into optimum magnetic hys- 
teresis damping constants, and the latter have to be obtained from computer 
solutions. The result of this optimization work makes possible a better 
design of the satellite with lighter altitude control weight, shorter rod, lengths, 
and smaller earth-pointing error than previously reported in articles in the 
Bell System Technical Journal. 

I. INTRODUCTION 

The dynamics analysis by Fletcher, Rongved, and Yu 1 has shown 
that a two-body satellite will achieve an earth-pointing motion from an 
initial tumbling as a result of energy dissipation in the hinge joint 
through the relative motion between the two bodies. For a practical 
application, we shall consider the earth-pointing body to be like a 
dumbbell and the auxiliary body like a sheet, the two being connected 
to each other through a hinge mechanism of universal-joint type to 
allow a two-degree-of-freedom relative motion. When the satellite is 
in an earth -pointing motion, the auxiliary body is parallel to the local 
horizontal in its unstable orientation, and the two axes of relative mo- 
tions are aligned with the roll and pitch direction (see Fig. 1). The over- 
all aspect of such a passive gravitational attitude control system has 
been studied by Paul, West, and Yu 2 employing the magnetic hysteresis 
damping mechanism. Certain designs of the satellite in terms of moments 
of inertia, spring constants, and damping constants have been given 
in both Refs. 1 and 2 for an altitude of 6000 nautical miles (nm). They 
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Fig. 1 — Two-body satellite configuration. 

are found workable, but by no means "optimum," which means the 
achievement of the best over-all damping performance in tumbling 
motion, large-angle motion, torque-free librations, and forced librations 
under environmental torques of various frequencies. 

The present work concerns a parameter optimization of the afore- 
mentioned two-body satellite employing two types of damping: the 
linear velocity of viscous type, and the nonlinear magnetic hysteresis 
type. With damping of the viscous type, analytical results can be ob- 
tained for the librational motions which may be used as a basis of 
parameter design in the case of magnetic hysteresis damping. The 
employment of magnetic hysteresis damping even makes the (linearized) 
equations of librational motions highly nonlinear, such that analytical 
treatment becomes intractable and results have to be obtained by 
numerical means. 

II. GENERAL EQUATIONS OF MOTION 

The equations of rotational motion which have been derived in Ref . 1 
will be repeated here. The coordinate systems, as indicated in Fig. 2, 
are denned as follows. Let O-XYZ be a nonrotating frame, with its 
origin at the geocenter 0, its Z-axis passing through the perigee of the 
orbit, and its F-axis parallel to the orbital angular momentum vector. 
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Let So-xyz be an earth-pointing frame, with its origin S at the center 
of mass of the composite satellite, with z axis parallel to OS„ (the local 
vertical) making an angle \p with OZ and with the y axis parallel to OY. 
The body coordinates of body 1, Si-xiyiZi , are defined along its principal 
axes, with adjusted 1 moments of inertia (h , h , h). Euler parameters 
(£> V, f> x) are employed to describe the motion of «S'i-£i2/iZi relative to 
S -xyz. The matrix of transformation from the latter to the former 
frame is given as 

2 - v- - f 2 + x 2 2(£„ + fx) 

I". J 2(fc-fx) -¥ + v 2 - F + x 

2(8- + «x) 2(-& + ijf) 



2(tf - nx) 
2(fc + if) 
-?-*• + ?■ + x 2 J 



(1) 



Among the Euler parameters the relation £ 2 + j? + f/ 2 + x 2 = 1 holds. 




Fig. 2 — Coordinate systems. 
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The coordinates S 2 -x 2 y 2 z 2 are denned along the principal axes of body 2, 
with moments of inertia (h,h, h). The rotation of body 2 relative to 
body 1 is specified by an angle, a, about the Xj axis and then an angle, 
/?, about the y 2 axis. The transformation matrix from Si-xtfjiZi to S 2 -x 2 y 2 z 2 
is 



[/>.;] = 



|~Q3 


SaS/3 


-CaS/3 





Ca 


oa 


_S/3 


-SaC/3 


CaC/3 



(2) 



where C and S are abbreviations of cosine and sine, respectively. 

The equations of motion of the two-body satellite with the hinge 
joint situated at the center of mass of the auxiliary body and at the 
earth-pointing Z\ axis of the satellite body are : 



IiUj = (h — h)(<»>2f»>3 — Gn 2 n 3 ) + T r \ + Tdi , 
I 2 6>2 = (I a — /i)(u> 3 coi — Gn 3 ni) + (T r2 + T d2 )Cct 

— 2' c Sq:, 

I3U3 = (h — h)(uiu 2 — Gnin 2 ) + (T r2 + T d2 )Sa 

+ T c Ca, 
/ 4 ci 4 = (h - /e)(wBW6 - Gn 6 n 6 ) - (T r i + T A )Cp 

+ W, 
-^5^6 = (-^6 — Ia)(o}6^a — Gntfu) — T r2 — Tai, 
ltd* = (h - / B )(«4«6 - Gn 4 n 6 ) - (Tri + T dl )S0 

- T r C/3. 



(3a) 
(3b) 

(3c) 

(3d) 
(3e) 
(3f) 



In the above, (ui,ui,mi) are the components of the total angular 
velocity of body 1 along S\-XiyiZ\ , and («< , m , co B ) are those of body 2 
along S 2 -x 2 y 2 z 2 . The coefficient G involves orbital elements: i.e., 



= 3ff(l - eT\l + eC+y, 

where e = orbital eccentricity, and 9. = 2t divided by the orbital pe- 
riod. Also, 



3 

I 



m = a l3 , tt l+3 = 2 &«*■ > * = !> 2 > 3, 
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with a,; and /;,-/ the elements of the transformation matrices (1) and 
(2). The constraint torque, T c , is given as 

T /S 2 a , C 2 a , S 2 /3 , C 2 ^-' 

c_= Wt h i* h 

^[(J 3 - /i)(«iwi - Gnm) + (T r2 + 7' rf2 )Ca] 

^2 



- 9? [(h - / 2 )(«i«2 - GWfe) + (T r2 + T M )Sa] 

Is 

1 \ 
+ ^ [(/ 4 - / B )(o)4« 6 - GWta) - (Tri + 7' dl )S/3] 



(4) 



+ a(a> 2 Ca + waSa) - j8(co 4 C/3 + w 6 S/3) > . 

The restoring torques T r \ and 7' r2 acting on body 1 along the Xi and y 2 
axes, respectively, arc linear with the relative angle of rotation: i.e., 
T r i = ha, Tri — h&, where h and // 2 are spring constants. The damping 
torques acting on body 1, T,n and T,, 2 , along the X\ and y 2 axes, re- 
spectively, are defined in the following. For viscous damping, Td\ = da, 
and T,i2 = C 2r 3, where C\ and C\ are viscous damping coefficients. For 
magnetic hysteresis damping, the torque is dependent on the history 
of motion and is defined in regions I and III of Fig. 3 for Ta-i as 

T,n = 2V+ T di tL -Y L - t (5) 

as long as | 7' d2 \ < T d2 , where /§, 7\/ 2 are constants, and /? , Tdi are 
the values of /J, 7^2 when /3 last changed sign. After | 7 1 </ 2 1 reaches 
fa-i. then T,n remains at fdi as long as /3 does not change sign, as repre- 
sented in regions II and IV of Fig. 3. The magnitudes of /3 and 7'</ 2 will 
be given in Section IV, where the minor loops will be described. Accord- 
ing to the major loop in Fig. 3, no energy dissipation will result if the 
amplitude of oscillation is less then p\ Tax is defined by replacing /3 
by a and subscript 2 by 1 in (5). 

It can be shown that the Euler parameters and the relative angles of 
rotation are related to the w,'s, i = 1, ■ • • , 6, as follows: 

k = KxXi " fX 2 + t?X 3 ), (6a) 

V - J(rXi + xX 2 - M, (6b) 
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L-#-J 



Fig. 3 — Magnetic hysteresis damping torque vs amplitude, Tj« — /3. 
t = |(-i?Xi+ £X 2 + xfc), 

x = K-^. - **■ - ft»), 

a = — oji + a> 4 C/3 + w«S/3, 
/3 = — co 2 Ca — co 3 Sa -f- <«>6 , 



where 



and 



X,- = Ui — a,i2Y) i — 1>2,3, 



J, = 0(1 - eVU + eC^) 1 . 



(6c) 
(6d) 
(6e) 
(6f) 

(7) 

(8) 



If the cu.'s in (3a) through (3f) are treated as dependent variables, 
then (3a) through (3f ) together with (6a) through (6f ) form a system 
of 12 first-order equations in the 12 unknowns £, rj, f, x, «, P and o>, , 
i = l, • • • , 6. This system of equations has been programmed on an 
IBM 7090 for numerical solutions for any given initial conditions and for 
given dimensionless parameters: satellite moment of inertia ratios 
I 4 /I s f i = 1,3, • • , 6, spring constants kt/I$?, damping constants 
d = d/I£l (viscous damping) or fdi/I&?, % = 1,2, &,P (hysteresis 
damping). Some numerical results have been given in Ref. 1, while 
more numerical results are summarized in Section IV for parameter 
optimization in the case of hysteresis damping. 
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III. PARAMETER OPTIMIZATION WITH VISCOUS DAMPING 

In the case of librational motion, the Si-XiyiZi and *S 2 -.t 2 ?/ 2 2 2 axes 
oscillate about the earth-pointing rotating frame S -xyz with infinitesi- 
mal angles £1 , ??i , ft and £2 , 172 , ft , respectively. Because of the uni- 
versal joint constraint, we have a = £ 2 — £1 , /3 = tj 2 — m , and f 2 = f 1 . 
The position vector of the hinge in the Si-Xiy^ coordinate system is 
— Lz 1 , while in the Sz-xtiy&z it is zero. Let us assume that the orbital 
eccentricity, «, is small compared with unity: i.e., e is assumed to be of 
the same order as the infinitesimal angles, £1 , 771 , f 1 , a, and /?. Hence, 
from (8) one obtains ^ = ft + 2cQCQf + 0(e 2 ), ^ = !M + 2eSQt + 0(e 2 ), 
and ^ = — 2ttfSttt + 0(e 2 ). Equations (3a) through (3f) can then be 
linearized in the case of viscous damping to two sets of equations in 
pitch and roll-yaw librations. Upon transformation of the independent 
variable from t to \[/ = Qt, these equations become: 

^r 1 - \- T d2 + Spim ~ k-S = 2eS^, (9a) 
ay/ !■> 

VL + T& + 7 r « + <3P* + X* 2 )/3 + 3pt* = 2eS^. (9b) 
(/i/'- ayr h 

TF ~ r r * + d - <^ IT + *&& - &" = °. ( 10a > 

ayr /1 fli^ 

rf 2 a f/ 2 £i ,1m , n „ -j dft 
T79 1 T7? "t" 7- -* rfl + U — </2 1 TT 
# 2 # 2 /4 # (10b) 

+ (4g 2 + M fci)a + 4gs£i = 0, 

g + (i - /. - /.)f, - (/. + A) § - A J - 0, doc) 

where 

7\/l = Clti -77 , jT</2 = C 2 ft -y- , 

«^ ay/ 

fa = /i-i//ifi 2 , £2 = fa/I^f, 

X = / 2 // 5 , M = /1//4 , 

Pi = (/1 - /i)/7« , p 2 - (/* - /i)/7i , 

gi = (/a - h)/h , q 2 = (h - h)/h , 

/1 = (/1 + /. - / t )/(/, + /«), and /, = (7 4 + 7, - /,)/ 
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Note that (10c) is obtained by adding up the linearized versions of 
(3c) and (3f). 

3.1 Free Librational Motions 

For the free pitch libration, i.e., with e = in (9a) and (9b), the 
stability criteria have been obtained in Ref. 1. To ensure that the 
earth-pointing frame, S„-xyz, be the equilibrium position, it has been 
found that the spring constant should be larger than the following critical 
value : 

«,.__ *BL. (11.) 

API + Vi 

Thus, the pitch spring constant & 2 is chosen as 

& = „&.= _ 32&PL, (llb) 

Xpi + P2 

where a > 1. Similarly, the roll critical spring constant is found to be 

&•=--*** (12a) 

Mgi + 92 
and the roll spring constant is chosen to be 

(,.«&._ _J*I« (12b) 

where b > 1. The characteristic equations in the complex variable s 
for the free pitch and roll-yaw librations are, respectively, of fourth order 
and sixth order: 



s 4 + C 2 (l + X)s 3 + [3(pi + p 2 ) - 3a(X + l)pW(Xpi + p 2 )]s 2 
+ 3C 2 (Xp, + jh) 8 + 9pip 2 (l - a) = 0, 

and 

•_l/Wu ^„5 , r U-X) + X(1 -g 2 ) 2 + MX(l - gi) 2 

S + Ci(l + m) s + k — ; — 7^ r? i \ 

L M(X + 1) - X(/zgi + «a) 



13) 



4g 2 [g 2 - m(& ~ Dgi] + 4giUgi - (fr - l)g 2 ] 

Mgi + 92 



_, " (1 + /x)U ~ X) + X(l - g 2 ) 2 + jzX(l - gi)(l - ft) 

1 m(X + 1) - X( M gi + ft) 



+ 4 
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+ 4fo + *) + <i - *)* M( '" 1 g ! )+ U (1 " g2 \ 1 s ° 

m(X + 1) -X(Mgi + 9 2 )J 

"(m - X)g 2 [g 2 - /*(& - DflJ (14) 

+ qilnqi - (b - DgJK/x - X) + X(l - g 2 ) 2 | 
(m£i + 92) {m(X + 1) - X(/xgi + & ) 1 

bfiXqiqiil — gi)(l — g 2 ) 



— 4(b - l)gig 2 



(m?i + 92) {m(X + 1) - X(/igi + g 2 ) } 

pgig^d - 32) 1 

L - 92(1 - gi)[g 2 - m(& - D?iJ 



- (1 - gi)/iX 



(m9i + g 2 ) {m(X + 1) - X(/xgi + g 2 ) } _ 



2 

8 



, .* (m - X)(m<?i + 92) , 16gig 2 (/J - X)(l - b) = 

1 M (X + 1) - \(mx + ff.) * m(X + 1) - X(/igi + g 2 ) 

In deriving (14) the satellite body has been assumed to be axisymme- 
tric. The optimization problem for the free librational motions is such 
that by varying the parameters in the coefficients of (13) and (14) one 
gets the largest possible negative real part of the complex roots for the 
most slowly damped mode of librations. These parameters are C 2 , 
a, X, pi , and p 2 for the pitch motion, and Ci , b, p } qi , and g 2 for the roll- 
yaw motion. The parameter X contained in the roll-yaw characteristic 
equation is determined from the pitch optimization and is not treated 
as a varying parameter in the roll-yaw optimization. From the optimized 
values of the six parameters X, /z, p t , p 2 , q\ , and g 2 , one can determine 
the optimized inertia ratios Ii/h , i = 1,3, • • , 6. The value of 7 2 is 
chosen from the magnitude of the gravitational torque at a given orbi- 
tal altitude such that, in the presence of various environmental disturb- 
ing torques, the forced libration amplitude will not be larger than a 
specified value. 

For a two-body satellite with geometric configuration of a dumbbell- 
sheet combination, we have 

Pi =1, V- = — 1> Qi = 1> and q 2 = —1; 
the characteristic equations are simplified to 

* 4 + (l-r-X)C 2S 3 + 3a^±i S 2 

x - 1 (15) 

+ 3(X - l)C 2 s + 9(a - 1) = 0, 
with parameters C 2 , a, and X, and 



(16) 

2 



58 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1965 

j-(l + ,)(»- A) + 4X + 4U _ 1} 1 fts 3 

|_ M + X J 

| r t( M -X)(l+M) + 4»A + 4>(t-l) , i(b , ? 1. 
|_ (/* - 1)(.M + X) J 

+ 4C, (m-XKm-1) s + 16(m-X)(6-1) = 

with parameters Ci , b, and /x. These two equations have been pro- 
grammed on a digital computer for computation of complex roots with 
a,b = 1.0 to 2.0, Ci , C 2 = to 7.0, X = 1.25 to 6.00, and » = 2.0 to 
24.0. The case with a,b < 1.0 will give rise to positive real parts of the 
roots and is of no interest to us. 

The root loci of (15) and (16) are plotted in Figs. 4 and 5, respec- 
tively, with fixed spring constants (a = 1.2, b = 1.2) and varying 
damping constants and inertia ratios. The roots are plotted only in the 
second quadrant of the complex plane, as they are complex conjugates 
with negative real parts. In the case of critical and overcritical damping, 
roots degenerate into the negative real axis. In the pitch case, with 
X = 2/[3 - 5 J ] « 2.61804, the two distinct modes will, at C 2 « 1.2805, 
coalesce into a single point, — n = — Re(s) = -1.16, on the negative 
real axis, corresponding to a l/e damping time of 0.137 orbit at large 
t or \p. This is the result given by Zajac. 3 For other values of X, there 
always exists a pair of complex conjugates for the two modes if d is 
not too large. It is noted from Fig. 4 that when X is in the range of 2.5 to 
4.0, a proper choice of the damping constant C 2 will make — Wi of the 
least damped mode not less than 0.7, which corresponds to a l/e damp- 
ing time of 0.23 orbit. Values of X in the above range and the corre- 
sponding optimum damping constants are given in Table I. This gives 
the satellite designer a wide choice of the inertia ratio, h/h , and the 
damping constant, Cz/I&t, and still the l/e damping time is not greater 
than 0.23 orbit in the free pitch libration case. 

From the complex root-locus plot for the roll-yaw free libration (Fig. 
5 at b = ki/ic* = 1.2), it is noted that the highest and the intermediate 
modes coalesce into a single point at n « 6.425, and Ci « 0.3625. 
The corresponding lowest mode has a poor damping. The intermediate 
and the lowest modes coalesce into a single point at n W 18 and Ci « 
0.0935, which gives a poor damping for the highest mode. A close exami- 
nation of the plot indicates that for n lying in the range of 8 to 10, 
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Table I — Optimum Damping Constant of Pitch Libration 



>(-£) 


C. 


Negative Real Part in $, — Re(j) 


Lower Mode 


Higher Mode 


2.5 
2.6 


1.425 
1.299 


0.88 
1.04 


critically 
damped 


2.61804 
2.75 


1.2805* 
1.165 


l.wt 

0.85 


1.16 
1.33 


3.0 
3.25 


0.995 
0.870 


0.81 
0.78 


1.17 
1.00 


3.5 
4.0 


0.80 
0.69 


0.74 
0.70f 


0.90 
0.75 



* Both modes coalesce into a single point on the negative real axis. 
t Re(s) = —1.16V 1 corresponds to 1/e settling time of 0.137 orbit. Re(s) = 
—0.70^ corresponds to 1/e settling time of 0.228 orbit. 

a proper choice of the damping constant, Ci , will make the negative 
real part of the least damped mode, — Re(s) m i n , larger than 0.24, 
which corresponds to a 1/e damping time of 0.66 orbit. If the lower 
bound of -Re(s) mi „ is relaxed to 0.20 (1/e time « 0.796 orbit), the 
range of n becomes much wider, i.e., 6 to 14, as indicated in Table II 
with the corresponding optimum damping constant. 

The variation of — Re(s) with the spring constants is such that the 
pitch optimum spring constant (for the least damped mode) depends 
on the choice of damping constant, C2 , at a given inertia ratio, X, as 
shown in Fig. 6 for a few chosen cases. In the roll-yaw case, the 

Table II — Optimum Damping Constant of Roll- Yaw Librations 



- (- '£) 


d 


Negative Real Part in ^-, — Re(s) 


Lowest Mode 


Intermediate 
Mode 


Highest Mode 


6 

7 
8 

9 
10 

12 
14 


0.65 

0.49 
0.375 

0.30 
0.25 

0.175 
0.13 


0.216 

0.233 
0.247 

0.258 
0.265 

0.280 
0.271 


0.211 

0.225 
0.250 

0.254 
0.250 

0.230 
0.212 


critically 
damped 

0.990 
0.715 

0.660 
0.495 



Note: Re(s) = —0.2V corresponds to 1/e settling time of 0.796 orbit. 
-0.28V corresponds to 1/e settling time of 0.568 orbit. 



Re(s) = 
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optimum spring constant for the lowest mode is found always smaller 
than that for the intermediate mode, as indicated in Fig. 7, while the 
highest mode gets better damping for increased k~i or b. For example, 
at n = 8.0, and C\ = 0.35, optimum k\ equals 1.13 k* for the lowest 
mode but equals 1.31 k* for the intermediate mode. From a practical 
consideration of the spring design, one should not choose the spring 
constant too close to the critical value because of possible decrease due 
to vibrations, thermal effects, etc. Furthermore, an analysis of pointing 
errors resulting from deviations in geometric configuration (due to 
rod deflections, 1 - 2 etc.) indicates the advantage of employing larger 
spring constants. In view of these conflicting results, one may have to 
settle for some compromise values, e.g., a,b tt 1.2 to 1.8. 

3.2 Forced Pitch Libration by Orbital Eccentricity 
The steady-state solution of (9a) and (9b) for 771 is found to be 

Vi = *W+ 6-VSV. (17) 

Here b\ = 2eA"'C,/ J 2 (/', - P->), and 
G, = 2eA-'[/VV' + (1 + 2\)k-J\P, 

+ &/V + (1 + X) (XPx + P,) {CU: 4- /c 2 2 )], 

where P\ = 3pj — 1, and P> = 3pz — 1. The ratio of the amplitude, 
jj, = (/<','-' + Oi), to the eccentricity, e, is plotted in Fig. 8 versus the 
spring constant, a, with varying €■• and X. This plot indicates the ad- 
vantage of using a lower spring constant, and, for 1.0 < a ^ 1.4, the 
advantage of employing a smaller damping constant. This latter ad- 
vantage is further reflected in the plot of rji/e versus C 2 (a = 1.2) 
(Fig. 9), especially when X is in the range of 2.5 to 3.5. It is noted from 
both Figs. 8 and 9 that at & = 0.5 (a = 1.2), rn/t is relatively inde- 
pendent of X. If the optimum damping constant for the free librational 
motion, as given in Table I, for X = 2.5 to 4.0 is used, the average value 
of rji/ 'e is approximately 2.25, which is only about 30 per cent larger 
than that given by C 2 = 0.5. 

IV. PARAMETER OPTIMIZATION WITH MAGNETIC HYSTERESIS DAMPING 

4.1 Energy-Fitting Method 

In the case of magnetic hysteresis damping, the equations of libra- 
tional motions are the same as (9a) and (9b) for pitch, and (10a) 
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0.8 1.0 



1.5 . 2.0 

a = k 2 /k 2 



2,5 



Fig. 8 — vi/* vs a with varying X and C> . 



through (10c) for roll and yaw, except that the damping torques Tai 
and Tdi become nonlinear and may be denned as follows. If, for example, 
the pitch relative angle of rotation from a neutral position, /?, is larger 
then j8, then the major hysteresis loop (see Fig. 3) will be traced, and 
(5) will be applied for the pitch damping torque. If however, /3 is less 
than /3 after a change of sign of /3, then the torque will in general trace 
a minor loop as, for example, that shown in Fig. 10. The maximum 
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Fig. 9 — rji/e vs C-> with varying X (a = k-</k->* = 1.2). 



magnitude of the torque in the minor loop may be a function of the angu- 
lar excursion. Since the torque always opposes the motion, the plot of 
the torque versus time is a broken curve (see Fig. 10). Each section 
of the curve can be represented by an analytical expression if it is ap- 
proximated as a section of a rhomboid or an ellipse. Such an approxima- 
tion can be used for numerical integration of the equations of motion. 

If the 7d2-/3 curve is approximated as a closed loop for a slowly 
damped system, the loop area which represents the energy of dissipa- 
tion per cycle of oscillation, En, is related to the angular amplitude of 
oscillation as follows: 



and 



E h = 47^(0 - 0) 



E h = KTtf* 



for 



^ 0o, 



for < 0o 



(18a) 



(18b) 



In the above, 0, O , K, and m are constants depending on the charac- 
teristics of the damping material, etc., where and /3 are found usually 
to be very small, e.g., = 1° to 4°, and 0o = 3° to 7°. If one equates 
En to the energy dissipation per cycle with the viscous damping, 
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E v = irCwp* (« = circular frequency of one of the principal modes) 
at a certain value of 0, then a relation between C 2 and f d2 results. For 
example, when E h = E v at /3 — ft, , then 



T d2 _ Vv /3o 2 u~\ 

htf " |_ 4 (ft> - /S) fij 



(19) 



A similar relation can be obtained between C\ and T dl /I^f. Thus, 
from the "optimum" viscous damping constants obtained in Section 
III, one may find the equivalent "optimum" hysteresis damping torque. 
Nevertheless, it is pointed out that the equivalent "optimum" T d 
obtained in such a way could be erroneous for the following reasons. 
First, this is not a slowly damped system, as can be observed from the 
computer solutions 1 of equations (3) with hysteresis damping, and also, as 
can be noted from Figs. 4 and 5, Re(s) is of the same order as Im(s) in the 
case of viscous damping. Hence, a closed-loop approximation for com- 
puting E h and E v based on a particular frequency is obviously a very 
poor one. Second, since the quadratic curve for E v can fit the i^-curve 
(m tt 1.5 power for /3 ^ ft and linear for ^ ft,) at only one point 
(see Fig. 11), (19) gives a much worse approximation at other values of 
(5 than at the point of fit. 

4.2 Numerical Method — Computer Solutions 

From the foregoing, it is apparent that "optimum" magnetic hystere- 
sis damping constant cannot be accurately evaluated from the "opti- 
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Fig. 11 — Ei,-P curve. 



mum" viscous damping constant by means of the energy method or 
other approximation methods, such as, the describing function method* 
(for example, see Ref. 4). Therefore, in the case of hysteresis damping, 
it is necessary to resort to numerical methods for the parameter optimi- 

* If one defines the equivalent gain as the ratio of the first harmonics of the 
output, hysteresis damping torque, to the amplitude of the oscillation angle, 
assumed to be sinusoidal, and equates it to CiS or C 2 s in the characteristic equa- 
tions, then the equations will only have terms of even power with, however, 
complex coefficients. This method is actually the same as the energy-fitting method 
and clearly offers no advantages. 
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zation. The results obtained in Section III for the viscous damping 
case may, however, be used as a guide. It is found from the computer 
solutions that the critical spring constants (11a) and (12a) obtained 
from the linearized equations (9a), (9b), and (10a) through (10c) 
with the viscous damping apply also for the case of hysteresis damping. 
It has also been demonstrated by computer solutions that the values 
of the inertia ratios, X and n, in the ranges given in Tables I and II, 
will give better damping for the same hysteresis damping torques and 
spring constants. We shall numerically integrate (3a) through (3f) 
and (6a) through (6f) with the hysteresis damping torque defined in 
(5) by putting a and /3 equal to 1° to 4°. This eliminates the unnecessary 
complex programming of the minor loops of the hysteresis damping, 
though no damping will result when a and are smaller than a and /3, 
respectively. If the initial condition is a tumbling motion, the numerical 
results will cover tumbling motion, large-angle motion, and librational 
motion. Only in the librational motion does the solution get less and 
less accurate when the relative angles get closer to 5 and /3. 

A number of computer runs have been made employing hysteresis 
damping constants in a wide range. The results indicate: (a) the number 
of orbits after which the satellite will stop tumbling from the initial 
condition, o>i = 0, ca = 50, m = 0, w 3 = m = o>6 = 0; and (b) the 
librational angles after 20 orbits. Computer runs have also been made 
(for the case of X = 3.333, n = 1A only) to determine the librational 
angles after 10 orbits from the initial condition: at t = 0, 6 X = 36°, 
6 U = 43°, 6 Z = 27°, a = 1°, = 10°, on = -0.840, <*> = 1.170, co 3 = 
-0.30. These angles are defined as X = C~\£ v &), dy = C \yi-y), 
and t = C^izi-z), (the caret denotes a unit vector), where d z is the 
earth-pointing error angle. The results of these computer runs are 
summarized in Figs. 12 through 14 for three sets of inertia ratios. 

It is noted from Figs. 12-14 that the damping constants which give 
better damping in librations are, in general, smaller than those for 
better tumbling damping. Poor damping will result in both tumbling 
and librational motions when f d i and f d i are either too small or too 
large. In the intermediate wide range of the damping constants there is 
relatively small variation in the damping time. The "optimum" damp- 
ing constants may be chosen from this wide range as, for example, 

0.19 - 0.38 for X = 2.62, M = 6.28 

W/ 2 2 = 0.16 - 0.32* for X = 3.33, n = 7.4 

0.08 - 0.16 for X = 4.0, y. = 14, 



* The upper limit is taken to be smaller than that given in Ref. 2. 
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and f dl is chosen to be one-half of f d2 . Since only a limited number of 
computer runs have been made, it is not possible to display T di /htt 
(i = 1,2) as a function of X(= 2.5 to 4.0) and M (= 6 to 14). The 
"optimum" damping constants corresponding to other values of X 
and n have yet to be determined from computer runs. 

The effects of spring constants have also been investigated on the 
computer by varying the spring constants under the same hysteresis 
damping constants. It is found that the variation in damping time is 
relatively small for ki/k* = 1.2 to 1.8 (t = 1,2). The lower limit is 
chosen to guarantee stability. 1 To ensure torsional fatigue strength 
of the torsion wire, it is preferable to choose high spring constants. 
Furthermore, employment of larger spring constants will reduce the 
error angle in case of rod bending, as indicated in the error analysis. 

The relation between the forced librational amplitude and the orbital 
eccentricity has been found from the computer runs to be 

2 £b 2e to 3c 

for e ^ 0.10. At c = 0.2, the satellite starts tumbling from an earth- 
pointing position after 1.5 orbits, whereas at t = 0.4 tumbling begins 
immediately after the start. 

Due to the nonlinear characteristics, the effectiveness of the magnetic 
hysteresis damping cannot be measured by the 1/e settling time as 
with viscous damping. However, an equivalent 1/e settling time for 
the hysteresis damping may be obtained in the following way. Take, 
for example, a typical computer run as plotted in Fig. 15 for 
2 [ X = 4, n = 14, T dl /htf = 0.06, T d2 /htf = 0.12, £< = lAk* (i = 1,2), 
a = p = 2°]. From the envelope curve of the amplitudes of the damped 
large-angle motion, the logarithmic decrements can be approximately 
evaluated by taking the motion as exponentially damped. It is found 
that the equivalent 1/e time ranges from one to three orbits with an 
average of two orbits. Within the validity of the librational motion, this 
represents an average 1/e time for both pitch and roll-yaw librations 
with the above parameters. 

V. SUMMARY — ILLUSTRATION OF A PRACTICAL DESIGN 

The results given in the previous sections indicate that there does not 
exist a single set of optimum parameters which give the best over-all 
damping performance. However, there have been found relatively 
wide optimum ranges of the inertia ratios, the spring constants, and the 
damping constants for the design of a two-body satellite. The optimum 
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NUMBER OF ORBITS 



Fig. 15 — Variation of 2 with magnetic hysteresis damping. 



ranges of the inertia ratios and the spring constants are found applicable 
to both viscous and magnetic hysteresis damping. They are X ( = I2/I5) 
= 2.5 to 4.0, n (= h/h) = 6 to 14, and ki, 2 /ici/ = 1.2 to 1.8. The 
optimum viscous damping constants, as tabulated in Tables I and II, 
can be read directly from the complex root-locus plots, whereas the 
optimum ranges of the magnetic hysteresis damping constants cor- 
responding to a few inertia ratios are obtained from the computer runs. 
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As far as the inertia ratios are concerned, it is advantageous to employ 
large X and n so that the deck rods 2 can be made shorter or the tip 
mass can be smaller. By employing shorter deck rods and furthermore 
by shortening the mast rod, 2 it is possible to balance out the solar torque, 
which was the largest disturbing torque for the satellite described in 
Ref. 2. As long as the solar torque is reduced to a small magnitude, the 
lowered gravity torque level resulting from the shortening of the mast 
rod will not substantially increase the error angles produced by other 
relatively small disturbing torques. In view of the foregoing, we could 
improve the design previously given in Ref. 2. The parameters may now 
be chosen as follows: 

I x = /, = 600, h = 10 slug-ft 2 

pitch gravity torque, T = 0.23 X 10~ 5 ft-lb/deg = 31 dyne-cm/deg 
X = 4, I 6 = 150 ;_/x = 14, h = 42.8 slug-ft 2 
fa/ltf = 0.06, T dl = 0.446 X 10" 4 ft-lb 
T d2 /m 2 = 0.12, T d2 = 0.892 X 10" 4 ft-lb 
k x = 1.8 k* = 0.55, A-, = 0.247 X 10~ 4 ft-lb/rad 
k 2 = 1.8 k 2 * = 1.8, ko = 0.804 X 10" 4 ft-lb/rad 
length of mast rod = 30 ft (assuming 160-lb satellite proper) 
length of deck rod = 16.5 ft (4 deck rods) 
roll tip mass = 2 X 2.5 = 5 lbs 
pitch tip mass = 2 X 9 = 18 lbs 

total mass for attitude control (excluding mast motor and support) 
= 26 lbs. 

The solar torque error angle is now only 1°, and the total error angle 
is found to be about 4° excluding the effect of the orbit eccentricity: 

magnetic dipole moment (TELSTAR's) 2° 

solar torque 1 

rod bending (silver-plated rod) 1° 

orbital eccentricity, e 2e to 3c 



total 4° + 2.5e. 

In the preceding, the error angle due to the residual magnetic dipole 
moment can be easily reduced to 1° or less by further refinement of the 
cancellation techniques used on the TELSTAR satellite. 
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